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Abstract 

A subset y C F" is a tile if F2 can be covered by disjoint translates of 
V. In other words, V is a tile if and only if there is a subset yl C F2 such 
that V + A — ¥2 uniquely (i.e., v + a — v' + a' implies that v — v' and 
a = a where v,v' &V and a, a £ A). In some problems in coding theory 
and hashing we are given a putative tile V , and wish to know whether 
or not it is a tile. In this paper we give two computational criteria for 
certifying that V is not a tile. The first involves impossibility of a bin- 
packing problem, and the second involves infeasibility of a linear program. 
We apply both criteria to a list of putative tiles given by Gordon, Miller, 
and Ostapenko in the context of hashing to find close matches, to show 
that none of them are, in fact, tiles. 



1 Tiles 

We first define tiles, and make some observations about them. Discrete tiles arise 
in problems in perfect codes jSj H] , and hashing [5] . Many of their properties 
have been extensively analyzed in [5]. 

In [5] the authors give a table of ten putative tiles, which is reproduced 
in Table [1] In Section 9 of [2] the authors give several criteria for non-tiling. 
However, none of these apply to the putative tiles given in Table [TJ Non-tilings 
have also been analyzed in [f and [6]. 

Definition 1. Let U he a finite- dimensional vector space over F2. A subset 
V C U is a tile (of U), if U can he written as the disjoint union of translates of 
V. In other words, there is a subset A (1 U such that every element of U can 
he written uniquely in the form v + a, where v & V and a G A. This may he 
expressed as 

V + A =U (1) 

{V + V)n{A + A) ={0}, (2) 

where X + Y -.^ {x + y : x £ X,y Y}. The equation ([TJ says that every 
element ofU can he written as a sum of an element ofV and an element of A. 
The equation ([2]) says that this representation is unique. Namely, v + a — v' + a' 
if and only if v + v' = a + a' ; thus uniqueness is equivalent to hoth sides heing 
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0. Note that the definition is symmetric in V and A, so that A is also a tile, 
called a complement ofV. 

Definition 2. If X C_ U denote by {X) the linear subspace generated by X (i.e. 
the smallest linear subspace containing X ). 

Note that is a tile of U if and only if F is a tile of {V). Namely, let W 
be a linear complement of (V) (i.e., is a linear subspace and W + (V) = U 
and W n {V) = {0}), and A be a complement of V in (F), then A + W \s a 
complement of V . If i? is a complement of V in U , then Br\{V) is a complement 
of V in (y). Thus we may assume, without loss of generality, that {V) = U. 
We say that is a proper tile of {V). 

If we also assume that the complementary tile A is proper, we say that the 
tiling is full-rank. 

Suppose that we are given a subset V C F'2 . We wish to know whether or 
not is a tile, and if it is, to find a complement A. In the next two sections 
we give two methods of showing that V is not a tile, and apply them to the list 
of putative tiles in Table [T] (taken from ^ using a slightly different notation), 
to show that none of them are, in fact, tiles. In this table, instead of listing 
bit vectors, as in [51, we list sets of indices, in which the corresponding bit 
vector must have a 1 . In the column labeled "generators of V" we give a list of 
generators. A set, S, is derivable from another set, T, if we can obtain it by one 
or more applications of either replacing an element with one of smaller index, 
or omitting it entirely. The set V will consist of the generators and any set 
derivable from them. Each of the putative tiles has cardinality 64. The column 
labeled k indicates that a putative complement A has cardinality 2*^. 



Table 1: Putative tiles 



k 


n 


generators of V 


6 


12 


{11}, {10, 5}, {9, 8} 


7 


13 


{12}, {10, 4}, {9, 8} 


8 


14 


{13,2},{13,1,0},{3,2,0} 


9 


15 


{14, 1,0}, {10, 2} 


16 


22 


{21,1} 


17 


23 


{22,0}, {19,1} 


18 


24 


{23,0}, {17,1} 


19 


25 


{24,0}, {15,1} 


20 


26 


{25,0}, {13,1} 


21 


27 


{26,0}, {11,1} 



For example, the first line in Table [T] lists generators 

{11},{10,5},{9,8}. 

Then V consists of the following 64 sets, each derivable from at least one of 
the generators: 
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{to}, < to < 11 : 

{10, to}, 0<to<5: 

{n, to}, < to < n < 9 
0: 



12 sets 
6 sets 
45 sets 
1 set 



2 Bin Packing 



A straightforward combinatorial approach to showing non-tihng is via bin pack- 
ing. This is the idea: Suppose that we're given a hnear projection 



For every w € W we have a "bin" 7r^^{w). We'h say that two vectors v,v' G V 
axe equivalent under tt (written v ~ v') if 7r(u) = n(v'). That is, equivalent 
vectors are always in the same bin. Since tt is linear, ii v ^ v' then v + a ^ v' + a. 
Thus we can lump all the vectors in an equivalence class together, and call this 
aggregate a "piece" . In order for V to tile U, it is necessary that we be able to 
pack together |J7|/|1^| copies of each piece so as to exactly fill up all the bins. 
More formally, we have 

Proposition 1. Let V Q U where U is a finite dimensional vector space over 
¥2, such that \V\ is a power of 2. Let n : U ^ W be a linear projection of vector 
spaces over ¥2, and c := \U\/\W\ (the bin size). For each integer i, < j < c, 
let hi be the number of equivalence classes of cardinality i, of vectors in V , 
where v ^ v' if and only if tt{v) — 7r(w'). Then a necessary condition for V 
to tile U is that there are non-negative integers dij (which are the number of 
pieces of size i that go into bin j ), such that 



Although a general bin-packing problem is NP-complete, quite often we can 
show non-solvability by simple arguments. 

This approach easily shows non-tiling for the last eight entries in Table [TJ 
In Table [2] below we give the results of projecting each of the putative tiles 
in Table [T] onto coordinates r, . . . , n — 1. The column labeled "piece census" 
specifies the multiset of piece sizes obtained. For example, in the row labeled 
fc = 8, the entry 10*5, 1*6, 1*8 means that there are 10 pieces of size 5, 1 of 
size 6, and 1 of size 8. We work out this example in detail to show the idea: 

There are ten pieces of size 5: 



tt:U ^W. 




(3) 



Y,d^,J ^ m/\V\)b^, for a^^^ = 0,...,c. 



(4) 



{{m, 2}, {to, 1, 0}, {to, 1}, {to, 0}, {m}} for to = 4, ... , 13. 



There is one piece of size 6: 



{{3, 2, 0}, {3, 1, 0}, {3, 2}, {3, 1}, {3, 0}, {3}}, 
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and one piece of size 8: 

{{2,1,0},{2,1},{2,0},{2},{1,0},{1},{0},{}}. 

For the rows with fc = 8 and k — 9, the bin size is 8, and the minimum 
piece size is > 4. Placing a piece of size 5 leaves no way to fill up a bin of size 
8. Similarly, in the rows for k = 16,. ..,21, the bin size is 4, and there are no 
pieces of size 1. Placing a piece of size 3 leaves no way of filling up the bin in 
which it is placed. 

Thus, none of the last eight rows in Tabled] can be a tile. 



Table 2: Projections 



k 


r 


bin size 


piece census 


8 


3 


8 


10*5, 1*6, 1*8 


9 


3 


8 


4*4, 8*5, 1*8 


16 


2 


4 


20*3, 1*4 


17 


2 


4 


3*2, 18*3, 1*4 


18 


2 


4 


6*2, 16*3, 1*4 


19 


2 


4 


9*2, 14*3, 1*4 


20 


2 


4 


12*2, 12*3, 1*4 


21 


2 


4 


15*2, 10*3, 1*4 



3 Linear Programming 

In this section we'll rewrite the defining conditions for a tile in terms of a 
linear program. We'll identify subsets S C U with their characteristic functions 
Xs '■ U — R: xs{x) — 1 ii X ^ S and otherwise. Denote convolution of 
functions /, g : i7 — ?• K by 

and the Fourier transform 

xeu 

As is well known: 

Note that if X, F C [/ then 

Xx^Xviz) := \{ix,y) : x e X,y e Y,x + y = z}\, 
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the number of ways of writing z as the sum of an element in X and an element 
of of Y . Thus, we may express the condition for F to be a tile (with A as a 
complement) as 

Xy*XA = l, (5) 

where 1 denotes the constant function with value 1. Although this is a necessary 
and sufficient condition (along with the condition that xa(u) G {0, 1}) it proves 
to be too weak to use as a linear programming criterion to certify non-tiling. 
We supplement it with the condition derived from ©: 

{xv * Xv)[XA ★ Xa) = \U\&, (6) 

where ^ : [/ ^ M is the function (5(0) = 1 and b{x) = when x 7^ 0. 
Taking the Fourier transform of ([5]) yields 

X^xl - \U\d. (7) 

This suggests a linear program in which we use the variables XA * Xa{u) 
instead of xa{u)- We are given V. Since = 2" we know \A\ (if it exists). 

We have variables 6„ for u S and Cx for x G . We'll want 



and 



We have the conditions 



bu = xa*Xa{u), 
Cx = \xa{x)\'^. 



< fou <|^| and is an integer, (8) 
< Cx <|^P and is the square of an integer, ([5^) 

60 =\A\ m 

Co dSb) 

6„ =0 if w ^ and xv * Xv{u) ^) 
Cx=Oiixy^O and xv{x) ^0 (gfe) 

Cx = ^(-1)'' "6„ for all x ([Sf) 

u 

If we drop the conditions about 6„ being an integer and Cx being the square 
of an integer, we get a linear program which must be feasible if y is a tile. 

One problem with this linear program is that it has a large number of nonzero 
coefficients. Just the condition that the Cx be the Fourier transform of the 
yields 2"(2" — max(| supp(xy ^Xv)!, I supp(xy)l) nonzero coefficients (we can 
either write the Cx as the transform of the bu or the as the inverse transform 
of the Cx, whichever yields a smaller system). We can immediately halve the 
number of nonzero coefficients by adding the equation 
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to the remaining equations for the Fourier tranform. However, the number 
of nonzeros is still quite large. We can greatly reduce this by means of ideas 
from the fast Fourier transform. We create new variables corresponding to the 
intermediate results of the transform. The usual sort of bookeeping now yields 
3n2" nonzeros (since each "butterfly" involves 3 variables) and introduces n2" 
new variables. 

Here are the details: Introduce variables t^.j with 0<i<n, 0<j< 2", 
with toj = Cj. For < i < n, < k < 2"", and < j < 2"^'^^, introduce the 
equations: 

^i+l,j+2"-»fe = iij+2"-*fc + ^ij"+2"-'fe+2i-»-i (9) 
ti+l,j+2"-^k+2"-^-^ = ti,j+2"-^k ^ j"+2"-'fe+2'»-'-i &) 

The values tnj are the values of the Fourier transform of Cj . We can also achieve 
a significant savings in our problem by noting that whenever one of the variables 
on the right hand sides of ([S]) is (which is the case for a significant fraction of 
the Cj), then we can "pass through" the remaining variable, or its negation, or a 
if both are 0, and not create a new variable. We note the effect of this special 
case in Table [31 by comparing n2" to the actual number of variables needed. 

One nice feature of the approach using linear programming is that the con- 
ditions like full-rank for the complementary tile A can be described as linear 
inequalities. 

Proposition 2. A subset A C U containing generates U as a linear subspace 
if and only if 

\XAix)\ < \A\-2, 

for all ^ X <E U . 

Proof. The value xa{x) is the sum of |^| terms each of which is ±1. Thus 
XAix) — \A\ if and only if x • a = for all a E A. This can happen if and only 
if A does not generate U. Similarly xa{x) — — \A\ if and only if x • a = 1 for all 
a € A. This is impossible since G A. Further note that xa{x) = |^| mod 2, 
thus the value of — 1 is impossible for |x^(a;)|. □ 

Thus, if we use variables representing |x(y4)p we can express full rank as 
|x(A)P<(|A|-2)2. 

For the first four of the ten putative tiles given in Table[T]the resulting linear 
programming problem was small enough so that either glpso30 or CPLEX[T] 
could handle it. The results of this approach applied to the first four cases in 
Tabic [His given in Table [1 

The system of linear equations for k ~ 8 were inconsistent. The linear 
programming problems for the last six rows of the table were too large for 
CPLEX to handle. 

^glpsol is the standalone solver contained in GLPK - the GNU Linear Programming Kit 
http : //www. gnu. org/sof tware/glpk 
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Table 3: Results from CPLEX 



k 


n 


n2" 


time in seconds 


rows 


variables 


nonzeros 


6 


12 


49152 


2.13 


33569 


33414 


99465 


7 


13 


106496 


214.33 


74349 


74710 


221693 


8 


14 


229376 


1.78 


140312 


142632 


419864 


9 


15 


491520 


269.71 


321016 


327828 


961832 



3.1 Farkas's Lemma 

Another nice feature of the approach using linear programming is the use of 
the Farkas lemma, which says that the infeasibility of a linear programming 
problem can be easily exhibited by means of a vector found from solving the 
dual problem. Thus, whenever the approach using linear programming is able 
to certify that a set is not a tile, one can find a vector, by solving the dual 
program given below, which is a certificate of non-tiling. 

Lemma 1 (Farkas). Let G,A,H and B be real matrices of dimensions mx n, 
s X n, m X r and s x r, and x,y, h and b be real column vectors of dimensions 
n,r,m and s. The system of of linear inequalities and equalities 

Gx + Hy>h (10) 
Ax + By^b (11) 

with a; > and y unconstrained, is infeasible if and only if there exist row vectors 
c > and d, of dimensions m and s respectively, such that 

cG + dA<0 (12) 
cH + dB^O (13) 
ch + db = l. (14) 

The two vectors, c and d produced by solving the dual program ()12|) consitute 
a certificate of infeasibility, which is easily checked. Calculating, we find that if 
PH)) is feasible then c{Gx + Hy) + d{Ax + By) > ch + db — 1. However, we also 
have c{Gx + Hy) + d{Ax + By) = {cG + dA)x < 0, which is a contradiction. 
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Abstract 

A subset y C F2 is a tile if F2 can be covered by disjoint translates of 
V. In other words, V is a tile if and only if there is a subset ^ C F2 such 
that V + A — ¥2 uniquely (i.e., v + a — v' + a' implies that v — v' and 
a = a' where v,v' £V and a, a' £ A). In some problems in coding theory 
and hashing we are given a putative tile V, and wish to know whether 
or not it is a tile. In this paper we give two computational criteria for 
certifying that V is not a tile. The first involves impossibility of a bin- 
packing problem, and the second involves infeasibility of a linear program. 
We apply both criteria to a list of putative tiles given by Gordon, Miller, 
and Ostapenko in the context of hashing to find close matches, to show 
that none of them are, in fact, tiles. 

1 Tiles 

We first define tiles, and make some observations about them. Discrete tiles arise 
in problems in perfect codes [31 13] , and hashing ^5 . Many of their properties 
have been extensively analyzed in [2]. 

In [5] the authors give a table of ten putative tiles, which is reproduced 
in Table [T] In Section 9 of [2| the authors give several criteria for non-tiling. 
However, none of these apply to the putative tiles given in Table [TJ Non-tilings 
have also been analyzed in [F and [7]. 

Definition 1. Let U be a finite- dimensional vector space over ¥2- A subset 
y C [/ (which we also refer to as a region^ is a tile (ofU), if U can be written 
as the disjoint union of translates ofV. In other words, there is a subset A C U 
such that every element of U can be written uniquely in the form v + a, where 
V G V and a G A. 

Note that the definition is symmetric in V and A, so that A is also a tile, 
called a complement ofV. 

Lemma 1. A subset V QU of a vector space over ¥2 is a tile if and only if 

V + A =U (1) 
{V + V)n{A + A) ={0}, (2) 

where X + Y := {x + y : x E X,y e Y} . 
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Proof. The equation ([T]) says that every element of U can be written as a sum 
of an element of V and an element of A. The equation ([2]) says that this 
representation is unique. Since we're working in characteristic 2, if v, v' £ 
V, a, a' € A then v + a — v' + a' if and only ii v + v' = a + a'. Uniqueness of the 
representation is equivalent to both sides being 0. □ 

Definition 2. If X <Z U denote by {X) the linear subspace generated by X (i.e. 
the smallest linear subspace containing X ). 

Lemma 2. A subset V QU is a tile of U if and only if V is a tile of {V) . 

Proof. Let be a linear complement of (V) (i.e., is a linear subspace and 
W + {V) a.TL<lW r\ ly) = {O}), and A be a complement of V in (F), then 
A + W \s a. complement of V . If _B is a complement of V in [/, then B n {V) is 
a complement oi V \a {V) . □ 

Thus we may assume, without loss of generality, that (y) — U. 

Definition 3. A subset V C U is a proper tile of U if V is a tile of U and 
{V) ^ U. 

If we also assume that the complementary tile A is proper, we say that the 
tiling is full-rank. 

Suppose that we are given a subset F C . We wish to know whether or 
not F is a tile, and if it is, to find a complement A. In the next two sections we 
give two methods of showing that V is not a tile, and apply them to the list of 
putative tiles in Table [1] (taken from f5| using a slightly different notation), to 
show that none of them are, in fact, tiles. 

In the following we freely use the common practice of identifying a bit vector 
with the set of indices in which the vector has the value 1. The motivating 
examples we use here are from calculations performed in [5] which is concerned 
with finding hash functions which are optimal in a particular sense. One way of 
constructing an optimal hash function H : ¥2 ^ ¥2 would be to find a region 
y C F2 with \ V\ — 2"^'" which is both a tile and optimal in the following sense: 

Definition 4. If S C F^\, and p e (0,1), define 

fsip) := E -P)r*^"^''^ 

a,beS 

where wt(x) is the Hamming weight of the bit vector x. Say that a region S C F'2 
is optimal at p if 

fs{p)>fs'{p) 

for all S' C F^, with \S'\ = 

We say that S is optimal if there is aQ < p < 1 such that S is optimal at p. 

Definition 5. A set system is a finite collection of finite subsets of the positive 
integers N. 
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We find it convenient to identify subsets S' C with the set system given 
by the non-zero coordinates of the vectors in S: 

S{S) := {{i e N : = 1} : w G S}. 

So in terms of set systems we have polynomials 

A,BeS 

where AAB denotes the symmetric difference. 

In order to state the theorem from [5] we define a partial order on subsets 
of N: 

Definition 6. We define a partial order on finite subsets o/N, called the right 
shifted partial order written S >r T , as follows: For a finite subset S* C N, 
denote by Si the i-th largest element of S. Then we say that S >r T if \S\ > \T\ 
and Si > Ti for 1 < i < \T\. 

Definition 7. We say that two set systems, S and T are isomorphic if there is a 
permutation n ofN which fixes all but finitely many elements, and a finite subset 
A G N, such that T = {AAt:{S) : S £ S}, where A denotes the symmetric 
difference. 

Note that if S is isomorphic to T then fs{p) — fT{p)- 

Theorem 1 (Gordon, Miller, Ostapenko). An optimal region S is isomorphic 
to an order ideal for the order >ji. 

In the table below we give the putative tiles taken from [5 . These were 
found by using the above theorem to restrict the search to all order ideals of 
cardinality 64. Instead of giving them as collections of bit vectors we specify 
them as set systems. Each is an order ideal for the order >jj which we specify by 
giving its set of maximal elements. Each of the putative tiles has cardinality 64. 
The column labeled k indicates that a putative complement A has cardinality 

For example, the first line in Table [1] lists generators 

{11},{10,5},{9,8}. 

Then V consists of the following 64 sets, each <_r at least one of the gener- 
ators: 

{m}, < m < 11 : 12 sets 
{10, m}, < TO < 5 : 6 sets 

{n, to}, 0<TO<n<9: 45 sets 
: 1 set 
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Table 1: Putative tiles 



k 


n 


generators of V 


6 


12 


{11},{10,5},{9,8} 


7 


13 


{12},{10,4},{9,8} 


8 


14 


{13,2},{13,1,0},{3,2,0} 


9 


15 


{14, 1,0}, {10, 2} 


16 


22 


{21,1} 


17 


23 


{22,0}, {19,1} 


18 


24 


{23,0}, {17,1} 


19 


25 


{24,0}, {15,1} 


20 


26 


{25,0}, {13,1} 


21 


27 


{26,0},{11,1} 



2 Bin Packing 

A straightforward combinatorial approach to showing non-tiling is via bin pack- 
ing. This is the idea: We are given a subset V C U of an F2 vector space and 
we wish to prove that it cannot be a tile. We find an auxilliary linear projection 

TT : [/ -> W^, 

where W is another F2 vector space. The projection tt partitions the elements 
of U into equivalence classes: w ~ u' if and only if n{v) = Tr{v'). We say that a 
subset of U of the form 

Pa,w {v + a : V eV, tt{v + a) — w}, 

where a £ U,w E W is a piece, and that a subset of the form Tr^^{w) is a bin. 
Note that Pa^w — Po.w+nia) + We then have 

Lemma 3. Let V d U be a tile with complement A <Z U , and n : U ^ W a 
linear projection of vector spaces. Then for all w £ W the collection {Pa.w ■ o, £ 
A,Pa,w ^ $} is a partition of the bin it~^{w). 

Proof. If Pa, HI n Pa', HI 7^ then there are v,v' £V such that a + v — a' + v'. By 
definition of a tile this implies that a = a',v = v' . We also have UaGAPa,™ = 
'!T~^(w) by surjectivity of tt □ 

This approach easily shows non-tiling for the last eight entries in Table [1] 
by using arguments about the cardinalities of the pieces Pa,ui- Note that since 
Pa,w = a+-Po,™+7r(a) that the possible cardinalities of Pa, u; are the same as those 
for Pq^u; In addition, by the above lemma, every one of the pieces Pq,™ ^ must 
be used in the partition of some Tr~^(w). In the arguments below we list the 
multiset of cardinalities of the non-empty Po,w for all w. 
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In Table m below we give the results of projecting each of the putative tiles 
in Table [T] onto coordinates r, ...,n — 1. The column labeled "piece census" 
specifies the multiset of piece sizes obtained. For example, in the row labeled 
fc = 8, the entry 10*5, 1*6, 1*8 means that there are 10 pieces of size 5, 1 of 
size 6, and 1 of size 8. We work out this example in detail to show the idea: 

There are ten pieces of size 5: 

{{m, 2}, {m, 1, 0}, {m, 1}, {m, 0}, {m}} for to = 4, ... , 13. 
There is one piece of size 6: 

{{3, 2, 0}, {3, 1, 0}, {3, 2}, {3, 1}, {3, 0}, {3}}, 
and one piece of size 8: 

{{2, 1,0}, {2, 1}, {2, 0}, {2}, {1, 0}, {1}, {0}, {}}. 

For the rows with fc = 8 and fc = 9, the bin size is 8, and the minimum 
piece size is > 4. Placing a piece of size 5 leaves no way to fill up a bin of size 
8. Similarly, in the rows for k = 16, ... , 21, the bin size is 4, and there are no 
pieces of size 1. Placing a piece of size 3 leaves no way of filling up the bin in 
which it is placed. 

Thus, none of the last eight rows in Tabled] can be a tile. 

Wc have not been able to find a projection to make this argument work for 
the first two rows of the table. 



Table 2: Projections 



fc 


r 


bin size 


piece census 


8 


3 


8 


10*5, 1*6, 1*8 


9 


3 


8 


4*4, 8*5, 1*8 


16 


2 


4 


20*3, 1*4 


17 


2 


4 


3*2, 18*3, 1*4 


18 


2 


4 


6*2, 16*3, 1*4 


19 


2 


4 


9*2, 14*3, 1*4 


20 


2 


4 


12*2, 12*3, 1*4 


21 


2 


4 


15*2, 10*3, 1*4 



3 Linear Programming 

In this section we'll rewrite the defining conditions for a tile in terms of a 
linear program. We'll identify subsets S CU with their characteristic functions 
XS '■ U R: xs{x) — 1 if X & S and otherwise. Denote convolution of 
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functions f,g:U-^M.hy 

f*9{x) = '^f{y)9{x + y), 
and the Fourier transform 

xeu 

As is well known: ^ ^ 

/ * .9 = / 5 ■ 

Note that if X, F C J7 then 

Xx*Xy{z) := \{{x,y) : x e X,y e Y,x + y ^ z}\, 

the number of ways of writing z as the sum of an element in X and an element 
of Y. Thus, we may express the condition for to be a tile (with A as a 
complement) as 

Xv*XA^l, (3) 
where 1 denotes the constant function with value 1. Although this is a necessary 
and sufficient condition (along with the condition that xa{u) G {0, 1}) it proves 
to be too weak to use as a linear programming criterion to certify non-tiling. 
We supplement it with the condition derived from ©: 

ixv * Xv)iXA * Xa) = \U\5, (4) 

where ^ : J7 ^> M is the function (5(0) — 1 and 6{x) — when x 0. 
Taking the Fourier transform of ([3]) yields 

x^£i = (5) 

This suggests a linear program in which we use the variables XA * Xa{u) 
instead of xa{u)- We are given V. Since [T^H^j = 2" we know \A\ (if it exists). 
We have variables bu for u G F2 and Cx for a; G F2 . We'll want 



and 



We have the conditions 



bu = xa*Xa(u), 
Cx = \xA{x)f. 



< &u <|^| and is an integer, (6) 

< Cx <|^|^ and is the square of an integer, ([6K) 

bo ^\A\ m 

Co m 

&„ =0 if w ^ and XV * Xviu) ^0 (01) 

=0 if X ^ and xv{x) ^0 (jSh) 

C:, =E(-1)^ "6„ for all 2; ([Bf) 
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If we drop the conditions about bu being an integer and Cx being the square 
of an integer, we get a hnear program which must be feasible if V is a tile. 

One problem with this linear program is that it has a large number of nonzero 
coefficients. Just the condition that the Cx be the Fourier transform of the bu 
yields 2"(2" — max(| supp(xy * Xv)\, I supp(xy)l) nonzero coefficients (we can 
either write the Cx as the transform of the or the b^ as the inverse transform 
of the Cx^ whichever yields a smaller system). We can immediately halve the 
number of nonzero coefficients by adding the equation 

u 

to the remaining equations for the Fourier tranform. However, the number 
of nonzeros is still quite large. We can greatly reduce this by means of ideas 
from the fast Fourier transform. We create new variables corresponding to the 
intermediate results of the transform. The usual sort of bookeeping now yields 
3n2" nonzeros (since each "butterfly" involves 3 variables) and introduces n2" 
new variables. 

Here are the details: Introduce variables Uj with < i < n,0 < j < 2", 
with to.j = Cj. For < i < n, < fc < 2'', and < j < 2"~'"^, introduce the 
equations: 

J+2"-»A;+2"-*-i = ^i,j+2"-*fc ~ j"+2"-»fc+2'»-*-i (EK) 

The values tnj are the values of the Fourier transform of Cj . We can also achieve 
a significant savings in our problem by noting that whenever one of the variables 
on the right hand sides of ([7]) is (which is the case for a significant fraction of 
the Cj), then we can "pass through" the remaining variable, or its negation, or a 
if both are 0, and not create a new variable. We note the effect of this special 
case in Table [31 by comparing n2" to the actual number of variables needed. 

One nice feature of the approach using linear programming is that the con- 
ditions like full-rank for the complementary tile A can be described as linear 
inequalities. 

Proposition 1. A subset A C U containing generates U as a linear subspace 
if and only if 

\xXix)\<\A\-2, 

for alio ^ X eU. 

Proof. The value xa{x) is the sum of |^| terms each of which is ±1. Thus 
Xa{x) = \A\ if and only if x ■ a — for all a E A. This can happen if and only 
if A does not generate U. Similarly xa{x) = — \A\ if and only if a; ■ a = 1 for all 
a € A. This is impossible since G A. Further note that xa{x) = \A\ mod 2, 
thus the value of \A\ — 1 is impossible for |xa(2;)|. □ 

Thus, if we use variables representing |xA(a;)|^ we can express full rank as 
\xX{x)\'<i\A\-2r. 
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For the first four of the ten putative tiles given in Tabic [T]the resulting linear 
programming problem was small enough so that either glpso30 or CPLEX[l] 
could handle it. The results of this approach applied to the first four cases in 
Table [His given in Tabled 



Table 3: Results from CPLEX 



k 


n 


n2" 


time in seconds 


rows 


variables 


nonzcros 


6 


12 


49152 


0.46 


33569 


33414 


99465 


7 


13 


106496 


72.60 


74349 


74710 


221693 


8 


14 


229376 


0.54 


140312 


142632 


419864 


9 


15 


491520 


39.72 


321016 


327828 


961832 



The system of linear equations for k ^ 8 were inconsistent. The linear 
programming problems for the last six rows of the table were too large for 
CPLEX to handle. 

3.1 Certificates of non-tiling from Linear Programming 

For the use of the linear programming criterion, we use the following well-known 
criterion to produce a certificate of non-tiling. 

Lemma 4 (Farkas). Let A £ Matmxn(R),& G K"- Then exactly one of the 
following statements is true: 

1. There is an y G M" such that Ay — b and y > 0. 

2. There is a z e M™ such that z > and z < 0. 

We explicity describe the matrix A and vector b for our problem. A solution 
to the second alternative above will consistute a certificate of non-tiling. Note 
that since the linear programming problem in second alternative is homoge- 
neous, in practice we replace it by 

There is a z G M™ such that A^z > and b'^z = — c for some 
convenient c > 0. 

For convenience we write f{x) — XA * XAix) The rows and columns of A 
are divided into regions. The first region (1) of rows is indexed by elements 
X G U. The next region of rows, which is indexed by the symbols f{x) for x € U 
corresponds to x G U for which f{x) has a fixed value - either |A| for a: = 
(region (2)), or for those x with xv{x) *Xv{x) ^ (region (4)). The next 
region is similar to the last, is indexed by symbols of the form f{x), for which 
/(x)^ has a fixed value - either |Ap (region (3)) for a? = or for those x for 
which xv{x) ^ (region (5)). 

^glpsol is the standalone solver contained in GLPK - the GNU Linear Programming Kit 
http : //www. gnu. org/sof tware/glpk 
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The columns are divided into two regions - the first indexed by the elements 
f{x) and the second by elements f{x). 



A-- 



/(^) 

T 


/(^) 

\ 


1 






1 


1 




\ — 


... 1 / 



(1) 



(2) 
(3) 
(4) 



(8) 



The vector 6 is everywhere, except for 6j(o) — 1^1' ^/(o) ~ 1^1^- 

We denote elements of C2 (the n-fold product of the group of order 2) by 
and elements of its dual, by x. The vector z giving the Farkas certificate 
is denoted by 2. Note that we've eliminated the upper bounds in equation ([S]). 

The part of the matrix denoted by T is the matrix of Fourier transform on 
C2 . We can use the same idea as in (O to make the resulting set of equations 
much more sparse at the expense of introducting auxilliary variables. 

In terms of equations these are: 

For }Fz < 0: 

^/(o)l^l +^7(0)1^1' <0- (9) 

As remarked above we set the right hand side of this to any convenient negative 
constant. Column /(O): 

+ ^/(o) ^ (10) 

Column f{x) where x ^ V + V: 

Y x{x)zs + Zfi^x) >0. (11) 

Column f{x) where x ^ V + V: 

> 0. (12) 

Column /(O): 

' ^0 + ^7(0) ^ 0- (13) 
Column /(x), where xvix) ~ 0: 

- Z£ + zj.j,. > 0. (14) 
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Column /(x), where xv{x) ^ 0: 

- Z£ > 0. (15) 

In practice, we found that we could require that inequalities ([T^ . (fH)) and (|15p 
be equalities, and we could still find a certificate of non-tiling. This effectively 
eliminates the parts of z indexed by x. 

For the first and third rows of Table [T] the linear program solver CLfH 
found particular easy to describe certificates. For the second and fourth rows 
the certificates did not have an apparent structure which allows a compact 
description. 

The certificates below are vectors which have 0.0 everywhere except in the 
coordinates specified in the tables below. The part above the line corresponds 
to indices of the form f{x) and the part below to indices of the form f{x). An 
expression like (c, I) means a segment of all coordinates starting at c of length I. 
An expression like (c, I, s) has a similar meaning except that there is a stride of 
s, meaning the arithmetic progression c, c + s, . . . ,c + ms where m is the largest 
integer such that c + ms < c + I. 






-1024.0 


320 


1024.0 


640 


1024.0 


(192,64) 


1.0 


(384, 64) 


1.0 


(576, 64) 


1.0 


(768, 64) 


1.0 


(1216, 64) 


1.0 


(1408, 64) 


1.0 


(1600, 64) 


1.0 


(1792, 64) 


1.0 


(2240, 64) 


1.0 


(2432, 64) 


1.0 


(2624, 64) 


1.0 


(2816, 64) 


1.0 


(3264, 64) 


1.0 


(3456, 64) 


1.0 


(3648, 64) 


1.0 


(3840, 64) 


1.0 



Table 4: Certificate of non-tileability for the first putative tile 

How widely applicable are these methods? We have tried the linear pro- 
gramming method on all of the examples of non-tiles given in [5], and in all 
cases it has found a certificate of non-tileability. In [B] (3 the authors show that 

■^The linear program solver included in the COIN-OR package http://www.coin-or.org. 
See [6] 
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-8192.0 


(0,16384,2) 


1.0 



Table 5: Certificate of non-tileability for the third putative tile 

there are no full-rank tilings in dimension 8 and 9 as a result of the execution 
of a very long running computer program. It would be interesting to see if the 
linear programming method could supply certficates of non-tileability for the all 
the cases examined. 
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